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Abstract. We derive the generalized Fokker-Planck equation associated with the Langevin equation (in 
the Ito sense) for an overdamped particle in an external potential driven by multiplicative noise with an 
arbitrary distribution of the increments of the noise generating process. We explicitly consider this equation 
for various specific types of noises, including Poisson white noise and Levy stable noise, and show that it 
reproduces all Fokker-Planck equations that are known for these noises. Exact analytical, time-dependent 
and stationary solutions of the generalized Fokker-Planck equation are derived and analyzed in detail for 
the cases of a linear, a quadratic, and a tailored potential. 

PACS. 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 05.10.Gg 
Stochastic analysis methods (Fokker-Planck, Langevin, etc.) - 02.50.-r Probability theory, stochastic pro- 
cesses, and statistics 



1 Introduction 

Introduced just 100 years ago [T], the Langevin equa- 
tion has become one of the most important and powerful 
tools for studying noise phenomena in systems coupled 
to a fluctuating environment [2]. The main advantage of 
this equation is that it provides a physically transparent 
and mathematically tractable description of the stochas- 
tic dynamics of such systems. The Langevin approach is 
especially effective if the noise that describes the action 
of the environment on the system can be represented as 
a time derivative, in the sense of generalized functions, 
of a stationary process with independent increments on 
non-overlapping intervals. In this case the solutions of the 
Langevin equation belong to the class of Markov processes 
whose properties are well known (see, e.g., Refs. [3,4,5, 
[6]). The stationary processes with independent increments 
and zero initial state constitute a class of Levy processes 
[7]. For brevity, we call the Levy process, whose derivative 
produces a given noise, the noise generating process. 

A large variety of physical, biological, financial and 
other systems is successfully described by the overdamped 
Langevin equation, i.e., the first-order stochastic differen- 
tial equation. One of the main statistical characteristics of 
these systems is the probability density of the solution of 
this Langevin equation. If the noise results from the noise 
generating process, then the solution possesses the Marko- 
vian property and the probability density satisfies a closed 
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equation. Sometimes this equation is called the differential 
Chapman-Kolmogorov equation [8] . On the other hand, it 
is referred to as the Fokker-Planck equation for some par- 
ticular cases. Specifically, a Gaussian distribution of the 
increments of the noise generating process corresponds to 
the ordinary Fokker-Planck equation jHOIB], and heavy- 
tailed stable distributions to the fractional Fokker-Planck 
equation [9lll ^ fTT ] fl2 p i3 | [T4] . In order to capture these im- 
portant cases, we call the closed equation for the probabil- 
ity density that corresponds to an arbitrary distribution 
of the increments the generalized Fokker-Planck equation. 
We have shown recently [TS] that for additive noise the 
generalized Fokker-Planck equation can be represented in 
a unified way through the characteristic function of the 
noise generating process. Since the form and coefficients 
of this equation depend fundamentally on the distribu- 
tion of the increments of the noise generating process, it 
provides a useful tool for studying the effects of different 
noises. 

The aim of this paper is twofold. The first is to de- 
rive the generalized Fokker-Planck equation that corre- 
sponds to the overdamped Langevin equation driven by 
multiplicative noise with an arbitrary distribution of the 
increments of the noise generating process. The second is 
to solve this equation for specific cases and, on this basis, 
to study the distinctive effects of different noises on the 
system. 

The paper is organized as follows. In Sect. [5] we discuss 
the overdamped Langevin equation driven by multiplica- 
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tive noise resulting from the noise generating process. In 
Sect. [3] we derive the generalized Fokker-Planck equation 
associated with this Langevin equation in terms of the 
transition probability density and characteristic function 
of the noise generating process. To confirm the validity of 
this equation, we consider in Sect.[4]various noises, includ- 
ing Poisson white and Levy stable noises, for which the 
corresponding generalized Fokker-Planck equations are 
known. Several exact solutions of the generalized Fokker- 
Planck equation are derived in Sect.O In Sect. [6] we sum- 
marize our results. 



2 Overdamped Langevin equation 

The temporal evolution of the relevant degrees of freedom 
of dynamical systems that interact with a fluctuating en- 
vironment is often described by the (dimensionless) over- 
damped Langevin equation 

i(t) = f(x(t),t)+g(x(t),t)at). (1) 

For different systems the variables in this equation have 
different meanings, but to be concrete we will call x(t) 
[x(Q) = 0] a particle coordinate, f(x,t) — —dU(x,t)/dx a 
force field, U (x, t) an external deterministic potential, £(f) 
a random force (noise) resulting from a fluctuating envi- 
ronment, and g(x,t) a multiplicative noise term. Under 
certain conditions (see, e.g., Refs. [4lll6j) the 'real' noise 
£(f) with finite correlation time can be approximated by 
an idealized noise that effectively captures all the essential 
noise effects and turns x(t) into a Markov process, with 
the result that many of these effects can be described an- 
alytically. 

Because of the singular character of the idealized noise, 
equation ([I]) has to be interpreted with care. The start- 
ing point relies on the fact [17] that this noise is the time 
derivative, in the sense of generalized functions, of the 
noise generating process T](t). According to this, the in- 
crement 5r/(t) = rj(t + t) — r)(t) of 77(f) is defined as the 

time integral, Srj(t) — L T dt'£(t'), in the sense of conver- 
gence in distribution. Therefore, the increment 5x(t) = 
x(t + t) — x(t) of the particle coordinate during a time 
interval r (r — > 0) can be written in the form 

Sx(t) = f(x{t), f)r + g(x(t), t)6r)(t), (2) 

which defines the meaning of equation (JTJ) in the Ito in- 
terpretation [18] (see also Refs. (20H1I6]). 

For a fixed r, the distribution of the increments 5rj(jT) 
(j = 0, 1, 2, . . .) is completely described by the transition 
probability density p(r]j+i,T\rjj), where 77^+1 and rjj de- 
note the possible values of r](jT + t) and r}{jr), respec- 
tively. These densities are properly normalized, /_ drjj+i 
x p(r]j+i,T\r]j) = 1, and satisfy the following condition 
lim T _^oP( ? 7j+i, T \ 7 lj) = S(Ar]), where S(-) stands for the 
Dirac S function and Ar\ = r]j + i — rjj. Moreover, if the 
first moment of 77 (jr + r) exists, it is assumed to be zero, 
i-e-, f™ oo dri j+l ri j+l p(r) j+1 ,T\r] j ) = 0. 



The noise generating process, i.e., the stationary Mar- 
kov process 77(f) = lim T _>o Ylj=o "V(J T ) w i tn ^(O) = 
and [f/r] denoting the integer part of f/r, is also com- 
pletely defined by the transition probability density 
p(r)j+i,T\r)j). We note in this regard that all transition 
probability densities of the form p(rjj + i 1 lT\i]j) (/ = 2, 3, . . .) 
can be expressed through p(r]j+i,T\r}j) by using the Chap- 
man-Kolmogorov equation [3H1[S]. In particular, for I = 2 
it yields p^^^r^) = dr] j+1 p(ri J+ 2,T\rj j+1 )p(T]j +1 , 
r\r)j). Thus, the statistical properties of solutions of the 
Langevin equation JT]) can be characterized by p(j]j + i , T\rjj) 
as well. Next, for simplicity, we additionally assume that 

If, for example, the transition probability density is 
Gaussian, i.e., 

e -A V 2 /(4Dr) 

^ T) -^mr^ (3) 

then 77(717-) = Y^j=o 3vU t ) ^ s a discrete-time Wiener pro- 
cess, which is fully characterized by two parameters, 

(Sri(jr)) = 0, (5r)(jr) S V (It)) = 2DS m t. (4) 

The angular brackets denote averaging over the incre- 
ments Stj^jt), and S nm is the Kronecker symbol. These for- 
mulas are the discrete-time versions of the mean (£ (f )) = 
and the correlation function (£(f)£(f')) = 2D5(t — t') of a 
Gaussian white noise of intensity D. Thus equations 
@ and (HJ) completely specify the Langevin equation (QJ 
driven by multiplicative Gaussian white noise [4]. 

To derive the Fokker-Planck equation, it is typically 
assumed that the first two moments of 6t)(Jt) exist [3j 
[5]. However, if p(Ar/,T) is a heavy-tailed function of Z\?7, 
then the second moment does not exist. In this case, the 
derivation of the Fokker-Planck equation that corresponds 
to the Langevin equation |T|) must be based solely on equa- 
tion ([2]). We emphasize that noises characterized by heavy- 
tailed transition probability densities p(Arj, r) differ qual- 
itatively from those characterized by p(Ari, t) with finite 
variances. Specifically, the latter have a frequency inde- 
pendent power spectral density dt e~ ZUJt {(;(0)t;(t)) (ui 
is the frequency) and consequently they are called white 
noises. In contrast, the power spectrum of the former does 
not exist. Nevertheless, they are a very useful tool for 
studying an important class of random processes that ex- 
hibit rare but large jumps. 



3 Generalized Fokker-Planck equation 

We define the probability density of the particle coordi- 
nate x(t) in the usual way: 

P(x,t) = (S(x-x(t))). (5) 

To derive the evolution equation for this probability den- 
sity, we need to be able to express the average values of 
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F(x(t)) and F(x(t), Sr](t)) in terms of P(x,t); the func- 
tions F(x) and F(x,y) are assumed to be determinis- 
tic. Keeping in mind the above definition of averaging, 
(F(x(t))) means averaging F(x(t)) over all increments 
Sr](jT) with j = 0, 1, . . . , [t/r] — 1 and r — > 0. It is ob- 
vious from equations ([2]) and ([5]) that the result can be 
represented as an average with respect to the distribution 
of x(t), i.e., 



with 



(F(x(t))) = / dxF(x)P(x,t). 



(0) 



In order to express (F(x(t),Si](t))) in terms of P(x,t), 
we use a two-stage averaging procedure [T5]. Since the 
variables x(t) and 5rj(t) are statistically independent and 
distributed according to the probability densities P(x, t) 
and p(Ar), t), respectively, we readily obtain 

/>oo />oo 

(F(x(t),Sr](t))) = dxP(x,t) dyF(x,y)p(y,T). 



To proceed, we introduce the Fourier transform, P k (t) , 
of P(x, t) according to the formula 



T{u(x)} = u k 



dxe %kx u(x), 



(8) 



and find, using definition (JSJ, that P k {t) = ( e - lfcr (*)). 
Ecjuation ((2]) implies that for t — > the increment of 
Pfc(t), i.e., <5Pfc = P&(f + r) — Pfc(t), can be written in 
the form 

5P fc = -tfcr(e- ite W/(z(t),*)) 

^_/ g -ifex(t)|- e — ikg(x(t),t)5n(t) _ ^g-j 

In accordance with ((6| and ((8]), the first term on the right- 
hand side of equation ([9]) reduces to 

ik(e- ik *Wf(x(t),t)) =F^f(x,t)P( X) t)}, (10) 

and equation ([7]) for the second term gives 

/ e -ifcx(i) f e -ikg(x(t),t)Sn(t) _ j\\ 



dj/e-^bw)^)-!]^^). (11) 



where p/c^) = P-tP^ T )} = (e~ lkSv ^) is the characteris- 
tic function of Sr](t). 

Substituting (jTUJ) and (jTT|) into equation ©, dividing 
it by r and taking the limit r — * 0, we obtain the following 
equation: 



^) = -W|/M)P(M) 



+ / d 2 /e- lfe V fcs (,,t)P(2/,i) (12) 



Una -\p k (r) 
t— >o r 



11 



(13) 



Since the transition probability density p(Z\?7, t) is nor- 
malized, i.e., Pq(t) — 1, the limit (|13p must satisfy the 
condition O = 0. If k ^ 0, then there exist three differ- 
ent cases, depending on how quickly p k (r) — 1 tends to 
zero as t — ► 0. First, if p k (r) — 1 = o(r), then 4> k = 
and the noise does not effect the system at all. Second, 
if Pfe(r) — 1 tends to zero slower than r, then \(f> k \ — oo, 
i.e., the influence of the noise is so strong that the sys- 
tem relaxes instantaneously to the final state. Finally, the 
case we are interested in corresponds to Pfc(r) — 1 = O(r), 
i.e., < \<f) K \ < oo and the noise acts on the system in a 
non-trivial way. 

We apply the inverse Fourier transform, defined as 



T 1 {u k } = u(x) 



1 

2^ 



dk e ikx u k 



(14) 



to equation (fT2| . Using the definition (fT3|) . we obtain 
where the function 



\g{y,t)\ \g{y,t) 



(15) 



(x) = lim — \p(x, r) 
t— »o r 



(16) 



is a special characteristic of p(Ar], r) for r — > that de- 
scribes the influence of noise on the system. Therefore, 
the desired generalized Fokker-Planck equation that cor- 
responds to the Langevin equation ([T]) driven by multi- 
plicative noise, which results from an arbitrary noise gen- 
erating process, takes the form 



^P(x,t) = -^-f(x,t)P(x,t) 



dy 



P(y,t) 
\9(y,WK9(y,t) 



(17) 



In accordance with the definition ([5]) , the solution of this 
equation must be normalized and satisfy the initial con- 
dition P(x,0) = S(x). 

To gain more insight into the connection between the 
generalized Fokker-Planck equation and the properties of 
the noise, we introduce the characteristic function S k = 
(e~ lfcl '( 1 - ) ) of the noise generating process rj(t) at t = 1. 

With the formula r)(l) — lim T ^o Sj=o ' 1 ^vU T )i it can 
be rewritten as Sk = lim r ^o(Pfc( T ))' 1 ^ T ' • Then replacing 
Pu{t) by 1 + T0fc and taking into account that lim e ^o(l + 
e) 1 ' 6 = e, we find Sk — e* fc , i.e., (j) k — \nS k . Thus, from 
equation (|12p we obtain an alternative representation of 
the generalized Fokker-Planck equation: 



d_ 

dt 



P(x,t) 



d_ 

dx 



f(x,t)P(x,t) 



~* 7 '{ F d y e ~ lkyp (y>^ lnS ^)}- 



(18) 
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In the particular case of additive noise, where g(x, t) = 1, 
equation (fT2|) becomes 



d_ 

at 



P k (t) = -F^f(x,t)P(x,t) } /',,( /W . 



(19) 



and the generalized Fokker-Planck equation l|18p simplifies 
to the equation 



-P(M) = - ^-/(a; ) t)P(a;,t)+^- 



l {P k (t)\nS k }, (20) 



which was derived in Ref. [15] . 

We note that the problem of deriving the generalized 
Fokker-Planck equation that corresponds to the Langevin 
equation |T]) has been considered earlier in terms of the 
Levy measure of the noise generating process [20,21,22, 
[23] . In contrast, the generalized Fokker-Planck equations 
([Tg|) and (|20[) are derived here in terms of the characteristic 
function S k of this process at t = 1. Since the stationary 
processes with independent increments are infinitely divis- 
ible, In Sk can be represented by the Levy-Khintchinc for- 
mula that connects Sk with the Levy measure [TJ. Hence, 
both approaches are, in fact, equivalent and lead to differ- 
ent forms of the generalized Fokker-Planck (see also Sect. 
14. 5[) . It seems, however, that the above approach which 
deals with the transition probability density of the noise 
generating process is more convenient for applications. 
Indeed, the transition probability density completely de- 
scribes the noise generating process and, in accordance 
with ([T3"]) and Sk = e^ fc , explicitly represents the charac- 
teristic function Sk- At the same time, there is no simple 
way to invert the Levy-Khintchine formula, i.e., to express 
the Levy measure through Sk [23] . 



4 Special cases of the generalized 
Fokker-Planck equation 

In order to confirm the validity of equations ([IT)) . (|18p 
and (|20p and to demonstrate their usefulness, we consider 
several specific noises for which the corresponding Fokker- 
Planck equations associated with the Langevin equation 
|T]) are already known. 



4.1 Gaussian white noise 

The transition probability density p(Arj, r) for Gaussian 
white noise is given by formula ([3]). Accordingly, Pfe(r) = 



-Drk 2 - 



-Dk 2 , and Sk = 



-Dk 2 



Then, taking into 



account that .F _1 -fPfc(*)* 2 } = -d 2 P(x,t)/dx 2 , we find 
that in the case of additive Gaussian white noise equa- 
tion ([20]) reduces indeed to the ordinary Fokker-Planck 
equation gU5J|B] 



d_ 

at 



P(x,t) 



d_ 

dx 



f(x,t)P(x,t) 



D^P(x,t). 



(21) 



If the Gaussian white noise is multiplicative, then 
-Dk 2 g 2 (y,t) and 



J kg(y,t 



dye lkv 4> kg {y,t)P{y,t) 



-DT- 1 {k 2 T{g 2 (x,t)P(x 1 t)}} 
D-^g 2 (x,t)P(x,t). 



(22) 



Applying the inverse Fourier transform to equation (|12p 
and using the above result, we again obtain the ordinary 
Fokker-Planck equation 



^P(x,t) = ~-^f(x,t)P(x,t) + D-^g 2 (x,t)P(x,t), 

(23) 

which corresponds to the Ito interpretation of the Lan- 
gevin equation ([T|) driven by multiplicative Gaussian white 
noise. We note that the last equation can also be derived 
from equation ([T7|) with (/>(■) — Dd 2 S(-)/dx 2 for this case. 



4.2 Poisson white noise 

As a second example we consider Poisson white noise, i.e., 
a random sequence of (5-pulses, defined as [251126] 



^(t) 



£(i) = J2*i8(t-ti). 



(24) 



Here n(t) is a Poisson counting process with the proba- 
bility P(n(t) — n) — (\t) n e~ xt /n\ of n > arrivals in 
the interval (0, t], X is the rate of the process, ti are the 
(random) arrival times of this process, and Zi are inde- 
pendent random variables of zero mean distributed with 
the same probability density q(z). It is assumed also that 
= if n(t) = 0. The noise generating process rj(t) 
is a step-wise constant Markov process whose increments 



Srj(t) = J f t+T dt'^it') are given by 



S V (t) 



0, ifn(r)=0, 
Xt T i^ ifn(r)>l. 



(25) 



In order to find the transition probability density 
p{Ar], t), we use the definition p(Ag, r) = (S(Arj — Sg(t))) 
which, together with |25|) . yields 



p(Ag, t) - P (r) 5(Ar)) + W{A V , r). 



(26) 



The first term on the right-hand side of this formula is the 
probability density of Sg(t) under the condition that none 
of the (5-pulses occurred during the time interval r. The 
second term, 

w(a v ,t) = y, p ^ ••■/ *M»»-£* 

n= l J —00 J — 00 \ „._1 / 



(27) 
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represents the probability density of Sr)(t) under the con- 
dition that at least one pulse occurs during this time in- 
terval. Taking the probabilities Pu{t) — P(n(r) — n) with 
linear accuracy in r, i.e., Po(t) = 1 — At, P\{t) = At, and 
Pu>2{t) = 0, we obtain from formulas (|2"6")l and (f2"7| 

p{Ar], t) = (1 - At) 6(Ar)) + \Tq(Ar]). (28) 

In accordance with the definition (|16p . for this proba- 
bility density 4>{x) — X[q(x) — S(x)], and the generalized 
Fokker-Planck equation (fl~7|) reads 

-P(x,t) = -—f( x , t )P(x,t) - \P(x,t) 



+A 



\g(v,t)\ \g(y,t) 



x-y 



(29) 



For g(x,t) = 1, i.e., in the case of additive Poisson white 
noise, this equation is consistent with those reported pre- 
viously [231I2%P?P8]. Of course, since S k = e^ 1 " 9 *), 
q(x) = J-~ 1 {qk} and 8(x) = .F _:L {1}, the same form of 
the equation (|2"9")1 (with g(x,t) = 1) follows from equation 
(|2"0|) as well. We note that a wide class of white noises, 
which are represented by random sequences of J-pulses 
with a mean number A of pulses per unit time, is char- 
acterized by the same transition probability density (j2"8"]) . 
Therefore, the generalized Fokker-Planck equation (|29p is 
also valid for all these noises. 



4.3 Compound noise 



Next we consider the noise £(i) = 5I m=1 £m(t) composed 
of a set of independent noises £ m (t) . In this case the noise 
generating process can be written in the form 

M [t/r]-l 

J7 (i) = lim V S Vm(jr). (30) 

m=l j=0 

Because of the statistical independence of the increments 
Sr] m (jT) of the partial generating processes rj m (t), the 
characteristic function S k = {e~ lkri ^) of ry(l) is expressed 
through the characteristic functions S mk = (e~ lk71m ^) of 

T] m (\) as follows: Sk = Ilm=i Smk- Therefore, in the case 
of additive compound noise the generalized Fokker-Planck 
equation (|2H)l becomes 

^- t P(x,t)=-^f(x,t)P( X> t) 

M 

+ J2F- 1 {P k (t)\nS mk }. (31) 

m— 1 

In particular, if M = 2 and £i(i) and are Gaussian 
and Poisson white noises, respectively, then equation (|31|) 
reduces to [5] 

poo 

-\P(x,t)+X dyP(y,t)q{x-y). (32) 



4.4 Levy stable noise 

The generalized central limit theorem implies that for 
a wide class of properly scaled transition probability den- 
sities p(Arj, r), the characteristic function S k corresponds 
to Levy stable distributions, S k — S k (a, (3, 7, p). It is well 
known [30) that S k {a, (3, 7, p) depends on four parame- 
ters: an index of stability a € (0, 2], a skewness parameter 
(3 £ [—1, 1], a scale parameter 7 G (0,oo), and a location 
parameter p G (—00,00). Assuming in accordance with 
the initial condition P(x,0) = S(x) that p = and ex- 
cluding from consideration the singular case when a = 1 
and 13 7^ simultaneously (in this case \4> k \ — 00), we 
obtain S k = <Sfc(a, /3, 7), where [30] 

S k (a, f3, 7) = cxp — 7|fc| Q (1 + i/3sgn(fc) tan — 

""(33) 

In the following we assume for simplicity that the con- 
dition g(y, t) > holds for all y and t. In this case 



(34) 



and the generalized Fokker-Planck equation (fT8"|) becomes 



d 



P(x,t): 



dt 

where 



d_ 

dx 



f(x, t)P(x, t)+T- 1 {G k (t) \nS k (a, (3, 7)}, 

(35) 



G k {t)=T{g a {x,t)P(x,t)}. (36) 

Equation (|35[) can be rewritten in a form containing 
the Riemann-Liouville derivatives defined as [31j 



s D° ± h{x) = 



(±1)" d n 
r(n — a) dx n 



dyh(xTy) y n ~ 



(37) 

where S D^_ and S D^ denote the operators of the left- and 
right-hand side derivatives of the order a (0 < a < 00), 
respectively. The function h(x) is defined on the interval 
[— s, s], n = 1 + [a], and r{z) is the Gamma function. 
Using the characteristic function (|33p . we first represent 
its natural logarithm as follows: 



In (a, /3, 7) = -7 



(l + /?)(ifc) a + (l-P)(-ik) a 
2cos(7ra/2) 



(38) 



Taking the Fourier transform of equation (|37[) with /i(x) = 
g a (x,t)P(x,t), we find 

(±ifc) a G fc (t) = ^{oo^g^x, i)P(z, *)}, (39) 

and combining this result with (|38|) we obtain the frac- 
tional Fokker-Planck equation 



^P(x,t) = ~f(x,t)P(,x,t)- 



7 



dt 



(1 + /3) 



dx 1 ' 2cos(7ra/2) 
3 £>" +(l-/3) OOJ D^].g Q ( 2 ;,i)P(x,t) 



(40) 



Equation (|40|) reproduces all known forms of the frac- 
tional Fokker-Planck equation that corresponds to the Lan- 
gevin equation ([!]) driven by Levy stable noise. It can be 
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easily rewritten in a form containing the Riesz derivative 
denned as I3T1 



d\x\ a 



h(x) = -T- 1 {\k\ a h k }. 



(41) 



With the help of this definition and the relations 



ira 



(oo^+ + oa D-)h(x) = 2cos ^- 1 {\k\ a h k } (42) 
and 

(oo£>? -ooD-)h(x) = 2sm^^- l {\k\ a - 1 h k }, (43) 

which follow directly from the Fourier representation 
ooD^.h(x) = T^ 1 {(±ik) a hk] of the Riemann-Liouville 
derivatives, equation (|40|) reduces to 



^-P(x,t) = ~f(x,t)P(x,t) + 7 -^-g a (x,t)P(x,t) 
at ox o\x\ a 



dx" 



+ 7/? tan — — — 



d 9 



(44) 



It is obvious that for a = 2 this equation takes the form 
of the ordinary Fokker-Planck equation (j23[) with D = 7. 
For a < 2, various special cases of equation (|44p have been 
given previously in Refs. [WHTOn^lT3llT4] . 



4.5 Infinite divisibility of the noise generating process 

It is well known (see, e.g., Refs. [71l2~4"]) that any stationary 
process with independent increments, including the noise 
generating process, is infinitely divisible. This means that 

the condition Sk — (S k )™ with being a characteris- 
tic function holds for each positive integer n. In this case, 
In Sk can be represented in the form [24] 



In Sk = / dz p(z)- 



-ikz 



1 + ik sin z 



(45) 



where p(z) is the density of the Levy measure of 77(1). Ac- 
cordingly, with the integral representation <$(■) = (l/27r) 
x dk e %k ^ of the S function, the last term in equation 
(fP5|) can be written as follows: 



dye- tk yp(y 7 t)\nS kg{y , t) 

dy[S[x-y- zg(y,t)} 



dz 



S(x-y) + sin z— g(x, t)S(x - y) J P(y, t). (46) 



Evaluating the integral over y and using the formula 
dyS[x - y - zg(y,t)]P(y,t) 

1 

d 



J2(- z r-g-^g n (x,t)p(x,t) 



n=0 



exp ( - z—g(x,t) )P{x,t), 



(47) 



which follows from the Taylor expansion S[x— y— zg{y, t)] = 
S^Lo( — l) n 9 n (y> t)(d/dx) n 5(x — y), we obtain the gener- 
alized Fokker-Planck equation 

§- t P(x,t) = ~f(x,t)P(x,t) 



dz 



Pi?) 



d 



exp I - z-^-g(x,t) 



-1 + sin z—g(x,t) 



P(x,t), (48) 
which was derived in Ref. |22j using a functional approach. 



5 Exact solutions of the generalized 
Fokker-Planck equation 

An important aspect of the generalized Fokker-Planck 
equation (|20|1 . which corresponds to the Langevin equa- 
tion (TjQ) driven by additive noise, is that in some cases it 
can be solved for all characteristic functions S k , i.e., for 
all noises represented by a time derivative of the noise 
generating process. This provides a unique opportunity to 
study in detail the effect of different noises on the same 
system. In Sects. 15.11 and I5.2[ we consider this problem 
for overdamped particles in linear and quadratic poten- 
tials, respectively. The stationary solution of the general- 
ized Fokker-Planck equation (fT5|) in the case of multiplica- 
tive Levy stable noise and a tailored potential is presented 
in Sect. lO 



5.1 Linear potential 

In this case U(x, t) — U(x) — — fox (/o is a constant force 
acting on a particle) and, since J r {df(x,t)P(x,t)/dx} = 
ikf Q Pk(t), equation (fT2"|) becomes 



^-Pk(t) = (-ikf + (t>k)Pk(t). 
ot 



(49) 



The solution of this equation, satisfying the initial con- 
dition -Pfc(O) = 1 that follows from the initial condition 
P(x,0) = 5{x), is given by 



Pfe(i) = exp(-ikf t + <t> k t). 



(50) 
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Therefore, using the relation S k = e^ k , we can write the 
solution of the generalized Fokker-Planck equation ([20)) as 



P{x,t)=F- l {Sl 



-iktf 



}• 



(51) 



In order to study how the behavior of free particles, 
when /o = 0, depends on the character of the noise, we 
calculate the second moment (x 2 (t)) = J_ tx dx x 2 P(x,t) 
of the particle coordinate x(t). Using the solution (|5ip 
and the integral formula S"(k) = — (l/2n) dx x 2 e lkx , 
we obtain 

(^(*))=-^ 5 L • ( 52 ) 
dk £ k=o 

Since in the case of free particles x(t) = rj(t), we have S k — 
-Pfe(l), and the normalization of P(x.t) yields So = 1- 
Moreover, if the first moment of x(t) exists, it is assumed 
to be zero, i.e., dSk/dk\k=o = 0. With these conditions, 
formula (l52l) reduces to 



(x 2 (t)) = (x 2 (l))t. 



(53) 



If the term (a; 2 (l)), which can be represented as 
(x 2 (l)) = lim T ^o<[Efio 1 " 1 ^0>)] 2 ) = lim r ^o<V(0))/r 
with (V(0)) - (Iv(t)-v(O)] 2 ) = IZo d yyMy,r), exists, 
i.e., (<5?7 2 (0)) = 0(t), then the noise is white and leads to 
normal diffusion of particles. In particular, (x 2 (t)) = 2Dt 
for Gaussian white noise, and (x 2 (t)) = Xt dyy 2 q(y) 
for Poisson white noise. If the transition probability den- 
sity p(At),t) has heavy tails, i.e., (5r] 2 (0)) — oo, then 
P(x,t) evolves in such a way that (x 2 (t)} does not exist. 
Such behavior is displayed, e.g., by the probability den- 
sity P(x,t) — J r ~ 1 {Sl(a, (3, 7)} of free particles driven 
by Levy stable noise with a £ (0,2) [9]. While correlated 
noise induces a variety of different diffusion regimes of 
free particles [32j . noise resulting from the noise generat- 
ing process can only give rise to normal diffusive behavior 
(if (8r] 2 (0)) < 00) or to non-diffusive behavior, character- 
ized by a probability density P(x, t) with infinite second 
moment (if (Sr] 2 (0)) = 00). 

5.2 Quadratic potential 

In the case of a quadratic potential, where U(x, t) = 
U(x) = bx 2 /2 (b > 0), equation ([HI) takes the form 

^Pk(t) + bk—P k (t) = P k (t)<f> k . (54) 
at ok 

Its general solution can be obtained via the method of 
characteristics [33], for example, and reads 



1 



P k (t)=exp\^-J dz^-+ Cl jV^-\n\k\+c 2 -tj, (55) 

where i'(x) is an arbitrary function, and c\ and c 2 are 
constants of integration. Since Pfe(0) = 1, the solution 
1(55)1 yields 

$(\\n\h\+o\ =expf-i / dz^-c\ (56) 



Replacing k by ke bt in this relation, we obtain 

-ke- bt 1 



^| fc | +C2 -t)=exp(-i^ 



dz Cj 

z 



rxy> 1 — — / (/:- 



Cl , 



and substituting this result into (|55p , we find that 



P k {t) = exp ( - ^ J dz- 



(57) 



(58) 



Finally, given that (j)k = hiSfe, the time-dependent solu- 
tion of equation (f20|) . P(x,t) = T^ 1 {P k {t)}, can be rep- 
resented in the form 



(59) 



We note that ]xi(S ze -bt /S z ) ~ —btzdhxS z /dz as 6 
0. Therefore 



1 



lim — / dz — In 



1 , S-n — b 



b J Z S z 



In Si 



(60) 



and in the case of free particles, i.e, 6 = 0, the solution ([59l) 
reduces to P(x, t) = T^ 1 {S t k ). This result is confirmed by 



the solution ([51)) with /o = 



5.2.1 Poisson white noise 



As a first application of the above results we consider Pois- 
son white noise. In this case <p k = — A(l ~ Ik), and formula 
(J5SJ) yields 



,-> , , / ^ f , 1ze- bt ~ Qz 

P k {t) = exp ( - - J dz 



(61) 



Next we assume that the probability density q(Arj) is ex- 
ponential, i.e., q(Arj) = (r/2) e -^l^l with r > 0. This 
implies that q k = r 2 (r 2 + k 2 )^ 1 , 



dz 



Qze- bt ~ Qz 



k r 2r l _ e -2bf )z 

dz 

(r 2 + z 2 )(r 2 + z 2 e~ 2bt ) 



1 , r 2 
= — In- 



2 r 2 + fc 2 e- 2b * ' 
and formula (161T) takes the form 



Pk(t) = 



r 2 + k 2 e- 2bt 



k 2 



(62) 



(63) 



For t — > 00 the time-dependent solution P{x, t) = 
J-~ 1 {P k (i)} of the generalized Fokker-Planck equation ([20)) 



8 S.I. Denisov, Werner Horsthemke, Peter Hanggi: Generalized Fokker-Planck equation: Derivation and exact solutions 



tends to the stationary solution P s t(x) = T 1 {P k (oo)}. 
Using the relation [34] 



dk 



cos kx 



(r 2 + k 2 ) s r{s)\2 



7T / \x\ 

r 



s-1/2 



K s -i/2(r\x\), 



(64) 

where Re s > and Ki (x) is the modified Bessel function 
of the third kind (or Macdonald function) [35] , we obtain 

/~2~ rlrlrl) 3 ^ 1 / 2 

P ^ x) = \l« 2>r{s) K -wWb (65) 

with s = A/ (26). This is the so-called i^-distribution, 
which is one of the basic distributions describing the sta- 
tistical properties of scattered waves [36]. It is interesting 
to note that for s < 1/2 the stationary probability den- 
sity exhibits singular behavior near the bottom (x = 0) of 
the potential well: P st (x) oc x ~ 2s if < s < 1/2, and 
P st (x) oc — In if s= 1/2. 



5.2.2 Levy stable noise 

In this case S k = S k (a, /3,7), and the characteristic func- 
tion ([33]) implies that 



In 



S. 



s z 



1 (1 - e- aM )\z\ a (l + z/3sgn(z)tan^) 



Using this result and the integral formula 



dz ■ 



sgnz 
1 



a 



sgn k 
1 



(66) 
(67) 



we find that 



K dz l\ n ^L = 2(i_ e -«»«)|fc|- 







S z a 

x ^1 + i/?sgn(fc) tan 
and the solution ([55]) reads 

1 - e- abt 
ab 



~2 



P(x,t)=T- L lS k [a,p, 7 



(68) 



(69) 



We note that special cases of this solution were known 
previously. Specifically, the time-dependent solution for 
symmetric Levy stable noise {(3 = 0) was derived in Ref. [S] , 
and the steady-state solution for asymmetric Levy stable 
noise in Ref. [371 . 



5.3 Tailored potential 

As a third example we derive the stationary probability 
density function P st (x) for overdamped particles interact- 
ing with a tailored potential 



U(x) 



c 1 dyyg a (y) 



(70) 



(c > 0) and driven by multiplicative Levy stable noise 
with g(x,t) = g{x) > 0. In this case f{x,t) — —cxg a {x), 
and equation ([35]) reduces to 



c—xG{x)+T- 1 {G k \nS k {a,l3, 1 )} = 0, (71) 

where G(x) = g a {x)P si (x). By applying the Fourier trans- 
form to equation (|71j) , we arrive at the ordinary differential 
equation 



ck— G k + G k In S k (a, /3, 7 ) = 0, (72) 



whose integration yields 



In— - = - [ dz~lnS z (a,(3,j). 
G c J Q z 



(73) 



Using the definition (|33|) and the integral formula ([67]) , we 
obtain 

G k = GoS k (a,f3,j/ac), (74) 



and so 



P st (x) =G g' a (x)^ 1 {S k (a,(3, 1 /ac)}. (75) 



To eliminate Go = (g a (x(t))) from the solution ([75]) . we 
use the normalization condition /_ dxP s t{x) — 1, which 
yields the desired probability density 



g- a (x)T- 1 {S k (a,P, 1 /ac)} 
JZ dxg-(x)^{S k (a,/3, 1 /ac)}' 



(76) 



To the best of our knowledge, the above result is the 
first stationary solution of the fractional Fokker-Planck 
equation l|44p. i.e., the generalized Fokker-Planck equation 
(jTBf associated with the Langevin equation ([1]) driven by 
multiplicative Levy stable noise. 



6 Conclusions 

We have derived a new form of the generalized Fokker- 
Planck equation associated with the Langevin equation 
for overdamped particles driven by multiplicative noise 
which results from the noise generating process whose in- 
dependent increments have an arbitrary distribution. The 
main advantage of this generalized Fokker-Planck equa- 
tion is that it accounts for the noise action in a unified 
way, namely through the characteristic function of the 
noise generating process at dimensionless time t = 1. Since 
the characteristic function is completely described by the 
transition probability density of the generating process, 
it is this density which ultimately determines the term in 
the generalized Fokker-Planck equation that describes the 
effect of the noise on the dynamics of the system. We have 
explicitly demonstrated this fact for various noises, includ- 
ing the Poisson white noise and the Levy stable noise. 

More importantly, we have solved the generalized Fok- 
ker-Planck equation in the cases of linear and quadratic 
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potentials driven by an arbitrary additive noise. A remark- 
able feature of these analytical solutions is that they give 
an opportunity to examine the effects of different noises 
on the same system. As an illustration, we have derived 
and analyzed in detail the time-dependent and stationary 
solutions that correspond to the Poisson white and Levy 
stable noises. Moreover, we have also presented the ana- 
lytical solution of the fractional Fokker-Planck equation 
that describes the stationary distribution of overdampcd 
particles in a specific potential driven by multiplicative 
Levy stable noise. 
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